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ABSTRACT 

Numerical simulations demonstrate that laminar breakdown in a boundary layer in- 
duced by the secondary instability of two-dimensional Tollmien-Schlichting waves to three- 
dimensional subharmonic disturbances need not take the conventional lambda vortex/high- 
shear layer path. 


1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the second author was in residence at the Institute for Computer Applications 
in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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The classical boundary-layer transition experiments of the early 1960’s [1], [2], [3] led 
to the perception that the laminar-turbulent transition process in a Blasius boundary layer 
may be viewed as a sequence of instabilities. The first two instabilities in the sequence, 
called the primary and the secondary instabilities, are linear at onset and are now very well 
understood. The primary instability in an incompressible flat-plate boundary layer is known 
as a Tollmien-Schlichting (TS) wave. The TS wave, as it equilibrates nonlinearly, becomes 
susceptible to what is known as the secondary instability. This is truly three-dimensional in 
character and is thoroughly reviewed by Herbert [4]. The known secondary instabilities are 
essentially of three kinds - the fundamental secondary instability, the subharmonic secondary 
instability, and the detuned secondary instability. In the first case, both the primary and 
the secondary wave have the same streamwise wavelength, and a streamwise vortex system 
on the scale of the primary wave (a “peak- valley splitting”) ensues in which the so-called 
lambda vortices are aligned in the streamwise direction. The nonlinear evolution of the 
fundamental secondary instability leads to detached high-shear layers. This stage suffers a 
tertiary instability manifested by a kink in the high-shear layer and an associated peak in the 
streamwise perturbation oscillogram, usually called a “spike”. Also associated with this stage 
are hairpin vortex elements. Thereafter, multiple spikes appear through what is possibly 
a rapid sequence of instabilities, and turbulent spot formation becomes imminent. This 
particular scenario, which is now known to have dominated all of the classical experiments, 
is called the K-type breakdown. 

The breakdown arising from a subharmonic secondary instability was observed exper- 
imentally much later [5]- [6]. In the subharmonic case, the streamwise wavelength of the 
secondary wave is double that of the primary wave, and the resulting pattern of lambda vor- 
tices is staggered. The detuned secondary instability involves so-called combination modes, 
i.e., a pair of secondary waves whose wavenumbers combine to yield the wavenumber of the 
primary wave. To date, high-shear layers and their associated spikes have not been observed 
experimentally for either the subharmonic or the detuned breakdowns. Instead, these in- 
stabilities appear to lead to breakdown by a rapid broadening of the spectrum, not just to 
the shorter wavelengths that are associated with the formation of “spikes” , but also to even 
longer streamwise wavelengths [7]. 

Numerical simulations of boundary-layer transition have produced the characteristic high- 
shear layers and spikes of the K-type breakdown [8], [9], [10], [11] and the early secondary 
instability stage of subharmonic transition [12]. In this paper we report numerical simulations 
which address the details of the laminar breakdown induced by subharmonic secondary 
instability. 

Numerical simulations of boundary-layer transition have been reviewed recently by Zang 
and Hussaini [9]. In the present simulations the parallel flow approximation is employed and 
the evolution is tracked in time rather than in space. The focus is on the vicinity of some 
point Xo and the approximation is that the displacement thickness is constant in x (although 
not necessarily in t ), the base flow is strictly in the streamwise direction and it is given 
by n b (y,t) = where u B (rj ) is the Blasius velocity profile which follows 

from the similar boundary-layer equations, and X — xo + c/t; Cf is the speed of a moving 
computational frame. Lengths are scaled by the displacement thickness, Sq, at Xo, and 
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velocities by the free-stream velocity, ttoo. In these units, X/xq = 1 + where d = 1.72. 

Although the boundary layer is assumed to be parallel, it is permitted to thicken in time 
[12], with the speed of the computational frame, c/, equal to the phase speed of the 2D wave. 
The numerical simulations reported here were performed with the boundary-layer version of 
the algorithm described in [13], with the nonlinear terms treated in skew-symmetric form 
[14]. 


If a and /? denote the fundamental wavenumbers in the streamwise and spanwise di- 
rections, respectively, then the fundamental wavelengths in these directions are given by 
L x = 2tt l a and L z = 2i r//3. The imposed periodicity lengths are 3 X L X and s z L z , where s x 
and s z are integers which specify the number of subharmonics that are permitted in each 
direction. (In the cases presented in this paper s z = 1 and 3 X = 1,2, or 9.) The ration- 
al numbers k x = k x js x and k z = k z /s z label the Fourier wavenumbers in the numerical 
representations with respect to the fundamental wavenumbers a and /3. 


2 


■mil 1 1 ml HI 


Case 

Re 

wave 

a 

P 

U) 

€ 

e 

*3 jp 

KKVF 

1100 

TS 

0.250 

0.000 

.08624 + ,00333i 

0.018 

0° 

l 


1100 

TS 

0.250 

0.209 

.09396 - .01263i 

0.008 



KKVS 

1100 

TS 

0.250 

0.000 

.08624 + ,00333i 

0.018 

0° 

2 


1100 

SQ 

0.125 

0.209 

.03385 - .01946i 

0.008 



KLSH 

730 

TS 

0.254 

0.000 

.09292 + .00160i 

0.00622 

36° 

2 


730 

SQ 

0.127 

0.241 

.03925 - .02254i 

1.78 x 10“ 6 



KLSL 

730 

TS 

0.254 

0.000 

.09292 + .00160i 

0.00622 

o 

CO 

2 


730 

SQ 

0.127 

0.241 

.03925 - .02254i 

4.5 x 10~ 6 



KLD 

730 

TS 

0.254 

0.000 

.09292 + .00160i 

0.00622 

36° 

9 


730 

SQ 

0.113 

0.241 

.03630 - ,02082i 

1.78 x lO" 5 




Table 1: Nonlinear Simulations 

The initial conditions for the simulations are of the form 

u(x,0) = &e{u fc (y,0) + e 2D e' e u 2D (y)e' ax (1) 

+j«3D u+ D (s,) + I £jd u - M 

where u 2 D(y) and u^jfy) are the least stable linear modes of the relevant equation (Orr- 
Sommerfeld or Squire) for the prescribed wavenumbers a and j3. These eigenfunctions are 
normalized so that their maximum streamwise amplitudes are 1, and, when used in the form 
given in Eq. (1), these maxima occur at x = 0 and z — 0. Thus, e 2 £> and e 3 £> measure the 
maxima of the streamwise fluctuations of the 2-D and oblique disturbances relative to the 
freestream velocity. The angle 6 measures the phase shift between the two-dimensional and 
the oblique waves. 

The parameters of five simulations are given in Table 1. They were chosen with particular 
experiments in mind. The table indicates whether the initial wave component is a TS 
wave or a Squire (SQ) wave. The temporal frequency, u t is included for reference. The 
direct numerical simulations were performed with sufficient resolution to meet the criterion 
recommended by Krist and Zang [15]; this typically required grids on the order of 192 x 96 x 
256 by the laminar breakdown stage. (However, spanwise symmetry was exploited to halve 
the number of grid points in that direction.) 

The first two simulations were performed under the conditions of the K-type experiment 
of Kovasznay, Komoda, and Vasudeva [2]. The KKVF and KKVS simulations addressed the 
fundamental and subharmonic laminar breakdowns, respectively, with the former simulation 
corresponding to the actual experiment. Figure 1 demonstrates that at these high initial 
amplitudes both types of breakdown exhibit strong detached high-shear layers. The figure 
shows the normal shear, du/dy, in the spanwise symmetry plane after 3| periods of the 2-D 
wave for the KKVF simulation, and after 4| periods for the KKVS simulation. In both 
cases a strong detached high-shear layer is present and by the time of the illustrations each 
has already undergone two roll-ups. (The figure corresponds to the two-spike stage.) Three- 
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dimensional flow-field visualizations of both cases reveal that the transition is associated 
with the formation of a large-scale lambda vortex/detached high-shear layer for which the 
breakdown occurs first in the symmetry plane and then spreads outwards in the spanwise 
direction. (The vortex lines for the KKVS case shown the left frame of Figure 2, for example, 
indicate the presence of a strong large-scale lambda vortex.) 

The subharmonic instability evolves more slowly in this case, as is typical of secondary 
instabilities at large amplitudes of the primary wave. Nevertheless, the KKVS simulation 
does furnish unmistakeable evidence that strong detached high-shear layers are not exclu- 
sively the province of K-type breakdowns. Continuation of these two simulations to later 
times shows that the strong similarity of the breakdowns remains, except, of course, that in 
the former case the patterns are aligned in the streamwise direction, whereas in the latter 
case they are staggered with twice the streamwise wavelength. We shall refer to this type 
of breakdown, which in this case is common to both the fundamental and subharmonic in- 
stabilities, as a primary-scale breakdown in deference to the scale of the vortical patterns 
(relative to the scale of the primary wave) which ultimately roll-up. 

The next two simulations were under the conditions of the Kachanov and Levchenko 
[6] experiments, which focused on subharmonic instabilities. The initial amplitudes here, 
particularly for the three-dimensional components, were much smaller than those for the 
other experiment. Indeed, the amplitudes were so low that laminar breakdown did not occur 
until after the primary wave had passed Branch II of the neutral curve. (In the KKVF and 
KKVS cases the breakdowns occurred well before Branch II was reached.) The KLSH and 
KLSL cases differed only in the amplitude of the initial 3-D perturbation. 

Figure 3 illustrates the time evolution of the normalized kinetic energy [9] in the 2- D and 
3-D components for these two simulations. On this scale Branch II is reached at t = 675, 
as evidenced by the subsequent decay of the primary wave. However, the primary wave 
eventually resumes growth because of feedback from the subharmonic waves [16]. 

Of more immediate interest for this paper is the character of the ultimate laminar break- 
down in these two cases. This is documented in Figure 4. The KLSL case does not break 
down in the same manner as the two simulations discussed previously. The symmetry plane 
contains only a very weak detached shear layer, and three-dimensional visualizations, such 
as the vortex lines for the KLSL case shown in the right frame of Figure 2, reveal that much 
smaller scale vortices form away from the symmetry plane. These are associated with some 
moderately strong detached shear layers. Eventually, strong shear does develop in the sym- 
metry plane, but it originates near the wall, as shown in Figure 5. Thus, this subharmonic 
breakdown does not follow the classical paradigm of a breakdown which commences in the 
symmetry, or “peak” plane and spreads outward in the spanwise direction. This breakdown 
is more diffuse, occurring at many spanwise locations, in some places in detached shear layers 
but in other places, such as the symmetry plane, in the near- wall region. We shall refer to 
it as a small-scale breakdown, since the vortices which eventually break down are on a much 
smaller scale than the primary wave. 

The KLSH simulation exhibits a stronger detached high-shear layer in the symmetry 
plane than the KLSL case. The shear layer of the latter case does eventually exhibit the 
familiar roll-up, although smaller scale features form off the symmetry plane as well. This is 
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an intermediate case between the distinct primary-scale breakdown of the KKVS simulation 
and the small-scale breakdown of the KLSL run. 

The principal factor determining the type of laminar breakdown which occurs appears to 
be the amplitude of the primary wave at the stage of incipient three-dimensionality, i.e., when 
the three-dimensional component reaches a level comparable to that of the primary. When 
the primary amplitude is large at incipient three-dimensionality, a primary-scale breakdown 
occurs; when it is small, a small-scale breakdown ensues; and when it is at an intermediate 
level, the breakdown has a mixture of both features. This speculation is supported by two 
additional numerical simulations. When the KLSL simulation is repeated, but with the 
boundary-layer thickness held fixed at its initial value, then the breakdown is manifestly 
a primary-scale one. Here the primary wave is always unstable and reaches a sufficient 
amplitude for a primary-scale breakdown to occur. Moreover, if an even larger amplitude is 
chosen for the KLSH case, for example, with e 4 = 1.78 x 10 -4 , then a clean primary-scale 
breakdown develops. Thus the small-scale breakdown of the KLSL simulation is associated 
with its small initial amplitudes and the eventual stability of the primary TS wave past 
Branch II. 

To date there have been no experimental observations of detached high-shear layers for 
subharmonic breakdowns. Such shear layers would be as easy to identify in the KKVS case 
as they were in the actual experiment [2] corresponding to the KKVF simulation. However, 
no experiments focusing on subharmonic behavior appear to have been done at such large 
amplitudes. 

At lower amplitudes experimental reports [6], [7] suggest breakdown by spectral broaden- 
ing to even longer scales and by various sum and difference interactions. The final simulation, 
KLD, allowed for both effects by leaving room in the computational domain for 9 stream- 
wise subharmonics of the primary wave and by setting the streamwise wavenumber of the 
initial three-dimensional component to be | the the streamwise wavenumber of the primary 
wave. Figure 6 shows the evolution of various harmonics. This simulation is indicative of 
a de-tuned, or combination, secondary instability which involves both the (4/9,1) and the 
(5/9, 1) modes. Note that although only the former is included in the initial conditions, the 
latter mode is quickly generated and both modes eventually evolve together. Notice also 
the rapid development of such other sum and difference harmonics as -y, and ^ . Flow 
field visualizations reveal that this particular simulation leads to a primary-scale breakdown. 
This type of simulation is much more expensive than those which allow for only one or two 
streamwise harmonics, so a more complete parameter study has not yet been performed. 
We suspect, however, that at sufficiently low amplitudes of the initial three-dimensional 
component a small-scale breakdown is also possible for the detuned instability. 

The simulations reported here aimed to illustrate that there are several paths to laminar 
breakdown in a boundary layer. The basic two types are primary-scale and small-scale, 
although features of both types can be present in a given case. 

The authors are pleased to acknowledge helpful discussions with Thomas Corke, Thorwald 
Herbert and Y. Kachanov. 
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Figure 3. Evolution of the kinetic energy in 
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